Applicability of Self-Consistent Mean Field Theory 
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Within the constrained Hartree-Fock (CHF) theory, an analytic condition is derived to estimate 
whether a concept of the self-consistent mean field is realized or not in level repulsive region. The 
derived condition states that an iterative calculation of CHF equation does not converge when 
the quantum fluctuations coming from two-body residual interaction and quadrupole deformation 
become larger than a single-particle energy difference between two avoided crossing orbits. By means 
of the numerical calculation, it is shown that the analytic condition works well for a realistic case. 
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I. 



INTRODUCTION 



The Hartree-Fock (HF) mean field theory has been ap- 
plied successfully to various quantum many-fermion sys- 
tems. By violating various conservation laws satisfied by 
the Hamiltonian, one obtains an appropriate mean field 
that could incorporate as many various correlations of 
the system as possible. Its many important theoretical 
concepts— such as a stability of the mean field and an ap- 
pearance of massless collective motion in association with 
restoration dynamics of broken symmetries— constitute 
an outstanding virtue of the self-consistent mean field 
theory. Here it should be noticed that the mean field is 
characterized by various averaged quantities, which ex- 
press an amount of symmetry breaking like a shape defor- 
mation in the coordinate space or a gauge deformation in 
the quasi-spin (pairing) space. These averaged quantities 
as well as the concept of mean field are expected to have 
physical reality only when the many-body system under 
consideration is near stationary, at least locally homoge- 
neous, and well isolated from the other local minimum. 
To the best of our knowledge, there has been no discus- 
sion of when a concept of the self-consistent mean field is 
realized or how and why it breaks down. In this paper, 
we will discuss to what extent the self-consistent one- 
body potential exists when two orbits around the Fermi 
surface interact with each other. 

Theoretically, an avoided crossing occurs rather fre- 
quently in various fields of physics like the molecular, 
atomic, biological systems, as well as the quantum dot 
and atomic nuclear systems 0, 0, El> 13 ■ I n a development 
of the nuclear structure physics, there have been many 
discussions on an applicability of the cranked mean field 
theory near an avoided level crossing region 0, 0, S @ ' 
where large angular momentum fluctuations and spuri- 
ous interactions between two crossing orbits have been 
explored. An argument for removing a certain spu- 
rious interaction and for introducing a set of diabatic 
single particle (sjd.) states in the cranked mean field 
HHIllIiimillEH seems to be reasonable, be- 



cause the angular momentum is a constant of motion and 
the interaction between different rotational bands should 
act at a given angular momentum rather than at a given 
rotational frequency. However, these treatments seem to 
be still in a phenomenological stage. 

The above argument may not be simply extended to 
deformation constrained mean field theory. In this case, 
an interaction between two potential energy surfaces 
(PESs) with different quadrupole deformation may not 
necessarily be regarded as spurious, because the Hamil- 
tonian and quadrupole operator do not commute with 
each other. Interestingly, it has been suggested that the 
pairing interaction between two orbits located below and 
above the Fermi surface and having spherical and de- 
formed shape, might be spurious in the level crossing 
region 9]. By eliminating the spurious interactions and 
comparing with experiments, very interesting conclusions 
have been deduced such as " a flat PES does not automat- 
ically lead to large fluctuation of the shape" 0. However, 
there still remains a decisive question of how the above 
statement is justified from the underlying dynamics. 

To understand how the mean field changes by itself 
when it acquires additional deformation, and to study 
more deeply what actually happens in the self-consistent 
mean field near the level crossing region from the un- 
derlying microscopic dynamics, a rather heavy numerical 
method called the configuration dictated CHF method 
[15J has been developed. Applying this method to real- 
istic cases, one may get many PESs which are approxi- 
mately characterized by the s.p. configurations relative 
to the lowest PES 0|. As pointed out in our previous 
papers 0, ^} , the CHF iterative calculation sometimes 
meets a difficulty of poor convergence or even noncon- 
vergence near the level repulsive region, no matter how 
much efforts one makes to get convergence. 

In this paper, we will explore how the competition be- 
tween the mean field and the quantum fluctuations com- 
ing from two-body residual interaction and quadrupole 
deformation plays a decisive role in breaking the concept 
of a self-consistent mean field near the level repulsive re- 
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gion. 



II. 



THEORETICAL METHOD 



The nuclear Hamiltonian is given by 

H = ^ Jaffna Cp + - ^] Vafosclt C^CgOy. (1) 



a/3 



Here and hereafter, a, /3, ■ ■ ■ are used to denote numerical 
basis. To define the self-consistent CHF state \^(q)) with 
a given quadrupole moment q, we apply the deformation 
CHF equation with quadratic constraints given by 



5((*(q)\H\V(q)) + ±w((*(q)\Q 2Q \*(q)) 



--a x (*(q)\£\*(q)) 2 )=0, 



(2) 



with 



<*(<z)|Q 2 ol*(<z)> - q, <*(«)|*|*(g)> = 0, (3) 



where x denotes a center of mass coordinate. Our nu- 
merical method of solving the CHF equation has been 
discussed in Ref. (HI, where the Gogny D1S interaction 
[II Gjl 113, IIJ 11113, the Coulomb force and the center 
of mass motion up to the exchange terms are taken into 
account. In Eq.Q, /i denotes an input parameter that 
allows us to vary an expectation value {^{q)\Q2o\^(q))- 
Meaning of parabola width w was discussed in Ref. [24J , 
and is chosen to be 1.0 x 10~ 3 McV/fm 4 in our calculation. 
The Lagrange multiplier A(g) is given by 



\(q)=w(fX-(V(q)\Q 2Q \y(q)}), 



(4) 



and an effective value of X(q) is allowed to change dur- 
ing the iterations. The symmetries Pe~ 17rJ ~ (z simplex) 
and Pe~ lvJy T 



(Sy) \2uL |26| are imposed in our numer- 
ical calculation, where P is the parity operator, e~' lnJi 
the rotation operator around i axis by an angle 7r, and f 
the time reversal operator. To keep the center of mass 
motion fixed, we impose a quadratic constraint in the 
cc-axis direction with Lagrange multiplier 



V = a x (y(q)\x\*(q)}, 



(5) 



and take a x — 1.0 x 10 _4 MeV/fm 2 in our numerical cal- 
culation. 

Having solved the CHF equation ©, one obtains a 
set of s.p. energies {ek(q)} as well as the s.p. states 
{(fik(q)}- Hereafter, the particle states are denoted by 
fx, v and the hole states by The letters k,l are 

used when no distinction is needed. To understand how 
the CHF state undergoes a structure change depending 
on the quadrupole deformation, it is desirable to obtain 
+ Aq)) in such a way that it can be regarded as a 



smooth function of q. For this aim, we apply the config- 
uration dictated CHF method, which is briefly recapitu- 
lated below. Let \^(q)) be a known CHF state satisfying 
condition (^(q)\Q2o\^(q)) = <?• To find a new CHF state 
\^(q+Aq)) that is supposed to be continuously connected 
with | ^ ((?)), we exploit the following condition 



lim (tpi(q)\ipj(q- 



A«z)> 



(6) 



where {fi{q)} denotes a set of occupied wave functions 
constructing the single Slater determinant \^(q)). That 
is, a small increment Aq is numerically adjusted by the 
maximum overlap criterion in Eq. © under a given ac- 
curacy, so as to maintain a characteristic property of the 
CHF state. In our calculation, Aq is so determined as to 
fulfill the condition 



\(<Pi(q)\<Pi(q + Aq))\ 2 > 0.9. 



(7) 



In this way, the configuration specifying \^(q)} is kept 
continuously as a function of q. Since the CHF state at 
q + Aq is dictated by the configuration of the preceding 
CHF state at q, this method is called the configuration 
dictated CHF method. It can also be generalized to get 
the excited HF states and the continuously-connected 
PESs. Applications of the configuration dictated con- 
strained Hartree-Fock-Bogoliubov (CHFB) method in 
the level crossing region as well as in the shape coex- 
istence phenomena, with and without the pairing effect 
have been reported elsewhere |16|. 

In our numerical calculation, a convergence condition 
is given in eV as 



E 



» 



(<?) 



>-i) 



(<?) 



< 10, 



(8) 



where e). denotes the s.p. energies in the nth iteration. 



III. 



FRAGILITY OF MEAN FIELD 



A. Difficulty of nonconvergence near the level 
repulsive region in CHF theory 

In our calculation, the s.p. wave functions are ex- 
panded on three-dimensional harmonic oscillator basis up 
to the principal quantum number Nq = 8. The ground 
state of 66 Se is obtained after having optimized triaxial 
deformation parameters of the Hcrmitc polynomials. The 
optimized range parameters thus obtained include some 
effects of higher major shells. Note that, no optimization 
has been done in the most HF and HFB calculations 
[27l I2H I29I l30l | , although a larger configuration space is 
adopted. It is well known that the optimum parame- 
ters change depending on the number of major shell, and 
become rather stable as the number increases [24j. In 
order to examine a reliability of the optimized configu- 
ration space, the ground state properties of 66 Se both in 
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TABLE I: The optimized ground state properties of 66 Se. 
The binding energy (BE) , quadrupole deformation parameter 
(P2) and triaxial deformation (7) are listed. The experimental 
BE is from Ref. 



(a) 





HF 


HFB 


Exp. 


BE [MeV] 


544.502 


545.623 


547.827 




0.241 


0.234 




7 (degree) 


47.256 


59.541 





HF and HFB calculations are listed in Table U together 
with the experimental binding energy |3lj . The binding 
energies in both calculations reproduce the experimental 
data well, though the total binding energy in the HFB 
is about 1.1 MeV lower than that in the HF calculation. 
A comparison between the HF and HFB calculations in 
Tab. [I] indicates that the nuclear deformation becomes 
small (favors the spherical shape) when one includes the 
pairing correlation into the mean field. In the following 
discussion, the number of major shell used in our calcu- 
lation may not have decisive importance. 

Starting from the ground state, the quadratic CHF cal- 
culation with the configuration dictated method is car- 
ried out in such small steps as the solutions are considered 
to be a continuous function of the deformation. Such a 
point-by-point heavy calculation is needed for discussing 
the dynamical structure change of nuclear system. In nu- 
merically obtaining the PES, the same range parameters 
as those in the ground state are used to trace an evolution 
of the ground state configuration as a function of defor- 
mation, which makes the s.p. level crossing dynamics 
transparent. Figure ^ shows (a) the quadrupole moment 
as a function of p, and (b) the lowest PES for 66 Se. One 
may observe that both the quadrupole deformation and 
PES change smoothly, except for a missing region around 
q = 150 fm 2 . When the constrained quadrupole moment 
is decreased by a small amount Aq from a critical point 
at qo = 177.365 fm 2 , the CHF iteration meets a diffi- 
culty of nonconvergence no matter how much efforts one 
makes to get convergence. After the missing region, the 
continuously-connected PES passing through an excited 
local minimum (q w 100 fm 2 ) is obtained. 

The nonlinear CHF equation is solved in an iterative 
way until the CHF Hamiltonian h{q) and density p(q) are 
diagonalized simultaneously such that [h(q), p(q)] = 0. 
In the numerical basis, the s.p. wave function tp^ 1 ' (q) at 
the nth iteration is expressed as <^afc an d P^(o) an d 
h( n \q) are given as 
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^v^pspffiq). 

■yd 



(9) 
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FIG. 1: CHF calculation for 66 Se: (a) the calculated 
quadrupole moment as a function of input quadrupole mo- 
ment parameter p; (b) binding energy as a function of 
quadrupole moment q. 



Here T(q) denotes the HF potential, and 

\™(q) = w(p - (^"-^l&ol^- 1 ^))), (10) 

where <F( n_1 ) (g)) is the Slater determinant constructed 

byl^iq)}. 

On a way to convergence, expectation values of 
one-body density p^ n '(q) in a representation of using 



M n) (?)},i-e 

Pkk 
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with respect to the hole and particle states are supposed 
to gradually reach 1 and 0, respectively. To explain what 
prevents the CHF iterative calculation from convergence, 
Fig. [21 depicts the diagonal components of proton and 
neutron densities as functions of number of iterations, 
for a case of nonconvergence with p, = 150 fm 2 (see 
Fig.^a)). In the case of nonconvergence, the quadrupole 
deformation parameter p is used in place of q. One may 
observe that the expectation values of proton density for 
the unoccupied orbits converge to 0, while those for occu- 
pied orbits to 1. For the case of neutron, there appears 
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FIG. 2: Diagonal components of (a) proton and (b) neutron 
density as a function of number of iteration for a nonconver- 
gent case with fi = 150 fm 2 . The single-particle basis where 
h is diagonal is used. Occ. and Unocc. stand for occupied 
and unoccupied orbits, respectively, (n, a) denotes the parity 
and signature, and its subscripts A and A + 1 represent the 
orbits responsible for the nonconvergent difficulty. 



a similar situation for the most single-hole and single- 
particle states, except for two specific orbits labeled as 
(—,—)a and ( — 7~ )a+i- Here, A denotes the number 
of neutron-occupied orbits, and the s.p. states are speci- 
fied by the parity and signature quantum numbers (it, a), 
since the asymptotic Nilsson quantum numbers are not 
good quantum numbers when the reflection symmetry is 
lost. For convenience, a subscript A and A + 1 for (tt, a) 
is introduced to identify the occupied and unoccupied or- 
bits responsible for the difficulty of nonconvergence. As 
seen from Fig. [3| these two specific neutron orbits with 
opposite quadrupole moments, i.e., one is of deformation 
driving and the other of anti-driving, lying just below and 
above the Fermi surface are interacting. One may expect 
that the two specific orbits (— , —)a and (— , — )a+i play 
a dominant role in preventing the CHF calculation from 
convergence. 

To exhibit crucial effects of two neutron orbits on the 
nonconvergence property in the CHF calculation, an ab- 
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FIG. 3: Neutron single-particle energies near the Fermi sur- 
face. Occ. and Unocc. stand for occupied and unoccupied 
orbits, respectively. The symbols are the same as those in 
Fig.H 



(«) 

A+1,A+1 



A, A 



and a difference of diagonal components h 

for the case with /i =150.0 fm are shown in Fig. 01 as a 
function of number of iterations. Here, the CHF Hamil- 
tonian is expressed in a representation where density ma- 
trix p( n > is diagonal, i.e., 



E(n- 
^afc 

a/3 



(q)h ( $(q)<P%- lJ (q). (12) 



>-i) 



(n) 



solute value of off-diagonal CHF Hamiltonian \h A 



A+ll 



In the CHF theory, it is usual to employ the above repre- 
sentation where the n-th quantities are entirely expressed 
in terms of the (n — l)th s.p. wave function 

In a case of convergence, the off-diagonal component of 
the CHF Hamiltonian would become smaller and finally 
reach 0, as the number of iteration increased. In the 
present nonconvergent case, there appears a staggering 
property in both the diagonal and off-diagonal compo- 
nents. During the iteration, (/i^ 1 ^, -J increases first and 
then starts to oscillate around some central value. In 
this case, the off-diagonal component always remains and 
never reaches 0, rather than being included into the mean 
field. Figure 0] indicates that the two interacting orbits 
make it difficult to apply the CHF mean field theory. 



B. Microscopic dynamics of nonconvergence 

This subsection discusses how the nonconvergent dif- 
ficulty and the staggering property appear as a result of 
the microscopic dynamics. Although the staggeringphe- 
nomena in Fig. 0] has been known as "ping-pong" [8| , it 
should be noted that the following analytic understand- 
ing and underlying physics are given for the first time. 
Namely, it is explored whether or not one gets a conver- 
gent CHF state at q = qo — Aq (a corresponding deforma- 
tion parameter /i = [1q — A/i), when there exists a conver- 
gent CHF state ^(qo)) at q = qo with corresponding set 



5 



0.30 

0.25 

.1 0.20 

c 

o 

I 0-15 
x 

g 0.10 
0.05 
0.00 



|hAA*ll 
hA+1A+1^hAA 




20 40 60 80 

Number of Iterations 



100 



FIG. 4: The absolute value of off-diagonal component and 
the difference of diagonal components of the CHF Hamilto- 
nian between two specific orbits as a function of number of 
iteration for the nonconvergent case with /j, =150.0 fm 2 . A 
representation of using s.p. basis where density p^ n ' is diago- 
nal is used. 



of self-consistent s.p. states {tkiqo), V?/c(<Zo)}- Since the 
CHF state |^(gn)} is used as an initial trial wave function 
to proceed the iterative calculation at q = qo — Aq, the 
following relations hold: 



(13) 



To study dynamical change of the s.p. wave functions 
during the iterations, it turns out to be preferable to 
use a fixed representation, i.e., the qo- representation of 
using the CHF s.p. states {ifik(qo)}, rather than the usual 
representation of using n-dependent s.p. wave functions 
{¥>L n_1 ^(<Zo — Aq)} as in Eq. ifT^l) . The nth Hamiltonian 
in the qo-representation is expressed as 

h M (Qo ~ Ag) = flk(qo)h^(q - Aq)(p i(q o ) 

a/3 



= tkl + rj™(<fo - Aq) - X (n \q - Aq)Q kl 

for n ^ 0, 
h$(qo - Aq) = hu(qo) for n = 0, 



(14) 



where the matrix elements of kinetic energy, deformation 
operator and the HF potential in the go-representation 
are defined as 



a/3 



t u = ^2<f* a k(qo)t a i3<fpi(qo), 

Qkl 

Aq) 



J2^k(qo)Q al 3<Pf3l(qo) 



a/3 



J2^k(qn) T ap(qn - &q)m(qo), 

a/3 

(15) 



and 

\W(q -Aq)=w{ fM) -An 

(q - Ag)|g 2 o|* ( "- 1) ('Zo - Ag)», 
for n^O, 

\W(qo - Aq) = X(q ) = w(n - <*(<Zo)|Q 2 o|*(9o))), 
for n = 0. (16) 

Making the following analytic understanding of our nu- 
merical results transparent and simple, we exploit such 
an approximate expression, i.e., 



A-l 



A-l 



i=l 



^^(q Q -Aq)^*(q -Aq), 



i=l 



which states that a contribution to the mean field from 
the lowest (A — 1) number of hole-states is independent 
of small deformation change Afi as well as number of 
iteration n. It turns out that Eq. i|17|) is well justified 
in our numerical calculation discussed in the previous 
subsection. Employing the above simplification, one may 
explore the nonconvergent dynamics governing the CHF 
iterative calculation in terms of a 2 x 2 truncated CHF 
Hamiltonian expressed as 



h^ A (Q0 Aq) 
/*a|i,a(<?o - Aq) 



h AA+l(lO - Aq) 
fcSi,A+x(«J - A?) 



(18) 



Although a set of eigen states {f k n \qo — Aq)} is numeri- 
cally obtained by diagonalizing the full CHF Hamiltonian 
^ (<7o — Ag), the characteristic feature of two interact- 
ing orbits is expected to be understood in terms of the 
truncated 2x2 Hamiltonian in Eq. (|18|) and 2x2 unitary 
matrix U^ n ' given by 



Aq) 



<Pa+i(Qo - Ac?) 



- t/ (n) 



9A{qo) 
<PA+i(qo) 



a (n) 5(«) 
c (n) 



(19) 



In accordance with using the go-representation, the above 
unitary transformation allows us to study the iterative 
process in terms of a mixing between two fixed states 
ipA^qo) and ifA+i(qo) irrespective of n. To explore a de- 
cisive role of relative phase between two orbits on the 
nonconvergent difficulty of the CHF iterative calculation, 
i.e., on the properties of resultant s.p. wave functions ob- 
tained after having diagonalizing h k ™\qo — Aq), we use 
four inter-dependent parameters in U^ n ' rather than a 
single independent parameter. Mixing parameters 

and dS n ' {br- n > = dP 1 ^ ) are a measure of the degree of 
mixing between the two specific orbits. In each diago- 
nalization, the subscripts A and A + 1 are used to assign 
the s.p. states in an energy-increasing order satisfying 
e^ ) (g -Ag)<e^ 1 (g -Ag). 
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With the aid of Eq. I|19f) and the approximate expres- multiplier \( n '> (qo — Aq) in Eq. I|16(l fulfills the following 
sion in Eq. I|17|) . it is easily shown that the Lagrange simple recurrence relation 



A(™ +1 > (q - Aq) = A< n > (q - Aq) + AA (n) (<fo - Aq) for n 0, 



where 



AX ( - n '>(qo~Aq)^~2wa^b^Q A , A+1 ~wb^{Q A+hA+1 -Q A ^} for n ^ 0, 

with an initial relation given as 

\W(q -Aq) = u;(/io-A/i-(*( )(go-AQ)|Q 2 o|* (0 )(go-Ag))) 
= X^iqo - Aq) - wAfi = \(q ) ~ wA[i for n = 0. 



(20) 



(21) 



(22) 



In the same way, one may derive the general expression of the matrix elements of 2 x 2 truncated CHF Hamiltonian 
in Eq. (|18|l at the (n + l)th iteration as 



*%a+i = h^ A+1 + a^b^Hv A+1AAA+1 +2wQ 2 A A+1 )+wb^ Qa,a+i(Qa+i,a+i - Qa,a), 

h Al 1] = h AA + 2wa (n '>b^Q AA Q AA+1 - {va+xaaa+x - wQa,a(Qa+x,a+x - Qa.a)}, 

^T+m+i = t$ +1>A+1 + 2wa [n) b {n] ] Q a +i,a+iQ a,a+i + & (n) {va+xaaa+x + wQ a +x,a+x(Qa+x,a+x 
for n^O, 

I 



M)}, 



(23) 



where the anti symmetrized two-body interaction in the 
go representation is defined as 

Vk x k 2 k 3 ki = ^2 ( Pak 1 {lo) l P* 1 k 2 { ( lo)v al pSfpk 3 {<lo)'PSk i (<lo)- 
afi-yS 

Here and hereafter, we use a simple notation h^ n+1 ^ in- 
stead of /i( n+1 ) (q a — Aq) etc, because we are only dis- 
cussing the CHF iterative process at q = qo — Aq. The 
matrix elements of truncated Hamiltonian hty are given 
as 

h A , A+1 = wAhQ a ,a+i, 

h { AA = e A (qo) +wAfiQ A ,A, 
h A+i,A+i = e A+i(qo) +wAfiQ A +i,A+i, (25) 



which are easily derived by using the initial condition in 
Eq. (JT3J) and initial relations in Eqs. lfH)l. l(T^|) and P^jl 
for the case with n = 0. In deriving Eq. I|25|l . the relation 



A, A n A,A+X 

l {0) h {0) 
l A+X.A n A+X,A+X 



£a(<7o) 







e A +i(qo) 



(26) 



which is satisfied by the preceding CHF solution at q — 
qo, is also used. 

With the aid of Eq. (|23|l , we get the following relations 
for two successive truncated CHF Hamiltonians ft! n ' and 

h (n+l). 



(n+l) 
l A,A+l 



(n+l) 
l A,A 



(n+l) 
l A+X,A+X 



= h { A ] A+ i + {a (n) 6 (n) - a^-W-VHvA+iAAA+i + 2wQ\a+i} 



+ w{b^ -b^ }Qa,a+i{Qa+x,a+i 

» 



Qa.j 



h<g A + 2w{a<ribW - a^-^ n -^}Q A , A Q A ,A+i 

{b (n) 2 _ b (n -^ 2 }{v A+1AAA+1 - wQa,a(Qa+1,A+1 

h A+i.A+i + 2itf{a< n W n > - o^-^-^Qa+i.a+iQa.ah 



!a,a 



)}• 



+ {b {n ^ 2 -b (n 1 ^}{v a +xaaa+x + wQa+x,a+x{Qa+x,a+x-Qa,a)}- 



(27) 
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In the same way, the off-diagonal component of the difference between and /i(™+ 2 ) is given as 



(n+2) 
l A,A+l 



An) 

L A,A+1 



+ { a (" +1 )&(" +1 ) - a {n -^b^}{v A , 



1AAA+1 



2wQ AA+1 } 



w 



{6« 



n+iy 



&(" ^ }Qa,a+i(Qa+i,a+i — Qa,a)- 

I 



(28) 
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FIG. 5: For the nonconvergent case with /i = 150 fm 2 (a) 
the off-diagonal Hamiltonian and difference of diagonal com- 
ponents as a function of number of iteration; (b) the relative 



phase a (n) 6 (n) of wave function for the orbit 



(n) 



Since the second term in the right-hand side of Eq. l(77|) 
contains a factor a^b^ whereas the third term has a 
factor frt™) 2 , the former term is retained in the following 
discussions because the parameter b^ is small. 

Numerical values of h^ A — h A L 1 A+1 and h^ AA+l are 
shown in Fig. HJa) for the nonconvergent case with fi = 
150 fm 2 . From this hgure, one may observe that both 
components exhibit a staggering property around some 
averaged values. Notice also that the off-diagonal compo- 
nent hj^A+i cnan g es its sign from iteration to iteration, 
whereas h AA — h A l_i A , j is always negative. The latter is 
easily understood because a couple of states tp A (q — Aq) 
and if A }_i(q — Aq) is defined so as to satisfy a relation 



6 a\io — Aq) < 6^^(170 — Aq). According to the per- 



» 



FIG. 6: For the convergent case with jj, — 180 fm . Notation 
is the same as Fig. |S] 



turbation theory of the 2x2 Hamiltonian, the sign of 
a (»5(™) for characterizing the lower state ip A \qo — Aq) 



(») 

A,A+ 



,(«) 

L A+1,A+1I 



111 



is the same as that of (h 

Fig. |S{b), the numerical value of oWftw is shown as a 
function of number of iterations. By comparing Fig.[SJa) 
with Fig. |SJb) , our numerical results are well understood 
in terms of the above simplified analytic expression in 
the two dimensional truncated space and the perturba- 
tion theory for the 2x2 Hamiltonian. Since the sign 
of a^b^ changes from iteration to iteration, it is clear 
that the higher state (p^^iqo — Aq) and the lower state 

(p A (qo — Aq) inter-change their properties from iteration 
to iteration. 

With the aid of coefficient a^b^ -a^" 1 )^ 1 ) in the 
second term of the right-hand side of Eq. Q27[l and from 
the above discussion on different signs between af^b^ 
and q^" -1 )^™ -1 ) ; one may understand why h^^ 1 ^ be- 
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comes large. Since there is a coefficient a^ n+1 'b^ n+1 ' — 
the second term of the right-hand side 
of Eq. (J2EJ, and since the signs of a ( Tl+1 ^ (n+1) and 
a (™- 1 )5(™- 1 ) are the same, it is easily recognized why 
the ph component at (n + 2)th iteration becomes small. 
These equations clearly state that the quantum fluctu- 
ations coming from two-body residual interaction and 
quadrupole deformation become small in one iteration 
and become large in the next, forming the staggering 
property. In other words, the major part of the two- 
body interaction could not be approximated successfully 
by the averaged one-body potential when the sign of 
a (")frO) changes from one iteration to the next. Physi- 
cally one may understand the above situation as follows: 
two mean fields, one characterized by occupied ipa(<1o) 
and the other by occupied <Pa+i (30)1 interact too strongly 
by the two-body residual interaction to be approximated 
by a single mean field. 

In contrast with the nonconvergent case, the same 
quantities for convergent CHF calculation at /i = 180 
fm 2 are shown in Fig. H3 One may observe that the sign 
of off-diagonal Hamiltonian stays the same during the 
iteration, and the sign of a^b^ is always the same dur- 
ing the iteration until the convergence is achieved. Since 
the sign of a^ n 'b^ n ' does not change during the itera- 
tions in the case of convergence, the quantum fluctuations 
characterized by the coefficient a^ n 'V- n ' — a^ n ~ 1 'b^ n ^ 1 ' in 
Eq. J27|) and by a (»+i)fe(«+i) _ a (n-l) h (n-i) in Eq J^J 

become small as n increases. That is, the quantum fluc- 
tuations are successfully incorporated into the mean field 
during the iterations. The above analytic and numeri- 
cal results indicate clearly that the sign of a(™)6(") be- 
tween two successive iterations is important in determin- 
ing whether the self-consistent CHF theory would be ap- 
plicable. 



IV. ANALYTIC CONDITIONS ON 
APPLICABILITY OF MEAN FIELD THEORY 

Since the analytic formulas within the two dimensional 
subspace discussed in the previous section nicely explain 



the characteristic feature of numerical iterative process of 
the CHF calculation near the level repulsive region, our 
next task is to determine what condition is capable of es- 
timating an existence of the mean field. In the previous 
section, it is clarified that the CHF theory breaks down 
when the off-diagonal Hamiltonian h < AA+1 changes its 
sign from iteration to iteration, whereas it is successfully 
applied when the sign is kept the same during the itera- 
tion. Namely, it is shown that the ratio of two successive 
off-diagonal components h^~Alil h^\ +l , [especially the 

first relation, h A ^ A+1 /h A V> A+1 ] plays a decisive role in the 
existence of the mean field. That is, the self-consistent 
mean field exists when 



■,(2) 

l A.A+l 
l A,A+l 



> 0. 



(29) 



From the discussion in the previous section, it is easily 
shown that the sign of h A n ~ A ^ 1 /h A n \ +1 is always posi- 
tive when the condition in Ea. (|29|l is satisfied, while it 
changes its sign when the condition is not satisfied. 

In this section, we focus on discussing what is meant 
by the condition 129|) . The ratio of two successive off- 
diagonal components satisfies the relation 



sign 



) n A,A+l 1 




j> = signj 


h (n) 1 
k ' l A,A+l J 





(n+1) 
l A,A+ 



3ign{/^V« H & (n) }, (30) 



because the signs of h AA+1 and a^b^ are always oppo- 
site, which is derived from our choice of states satisfying 



» 



(?o 



Aq) < e<r>(g 



Aq) and from the discussion 



based on the perturbation theory and is justified by our 
numerical calculation. With the aid of Ea. (|23[) . the last 
quantity in l|3U|) is expressed as 



fejgt+i = h% A+1 + a^b^(vA+iAAA+i + 2wQ% A+1 ) + wb^ 2 Q A ,A+i(QA+i.A+i - Qa,a) , 



Since an applicability of the mean field is supposed relation by using the following relation 
to be decided by a competition between the quantum 

fluctuation {ph component of hS 1 ') and the mean field , / £ ( 1 ) q \ 

part (pp and hh components of h^), let us find their U^h^U^ = J A ^ , (32) 



which should approximately hold within the truncated 
space under the assumption in Eq. Q1TJI. From an off- 
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diagonal component of the above relation, we obtain 

-h (1) -Z^ + ^W 1 ' (331 

'M+l.A+l A, A — \j(X) ^ (x) ] n A,A+V y° ) 
I 



Since the factor (M 1 ' 2 — a^ 2 ) is always negative [because 
the relation | a.^ 1 -' | > \b^ \ holds in the unitary transfor- 
mation in l|19|)l. one finally gets a relation, 

f h {2) } ( h {2) } 

a A.A+l n A,A+l I 

= sign{e A +i((7o) - £a(<Zo) - {^a+iaaa+i + 2wQ AA+1 ) 
—wAfi(Q A>A — Qa+x,a+x) — 0(b)}. (37) 

With the aid of this relation, our condition l|29l) is ex- 
pressed as 

eA+l(<?o) ~ £A(qo) > VA+1AAA+1 + 2wQ AA+1 

+ wAn(Q A . A - Qa+i.a+i) + 0(b), 

(38) 

which guarantees an existence of the self-consistent mean 
field, whereas the opposite condition, 

£A+l(go) - £A(qo) < VA+XAAA+X + 2wQ%A+X 

+ wA[i(Q A , A - Q A+ i <A+ x) + 0{b), 

(39) 

states a breakdown of the mean field. Note that the 
sign of the two-body interaction is important in Eas. (|38|) 
and l|39l) . which is in mark contrast with the well-known 
stability condition of the mean field. Furthermore, the 
present condition on a breakdown of the mean field has 
single-particle character, whereas the stability condition 
of the mean field has collective character. 

The physical meaning of the condition in Eq. I|38|l is 
clear: two-body correlation between nucleons can be suc- 
cessfully incorporated into the mean field as much as 



Inserting Eq. (|25|l for (|33|l produces the relation 



(34) 



(35) 



(36) 



possible when the energy difference between two inter- 
acting orbits is not smaller than the quantum fluctua- 
tions coming from the two-body residual interaction and 
quadrupole deformation. 

Fig. EJa) gives the s.p. energy difference between two 
specific orbits <p A (q) and <pA+i(q), and the quantum fluc- 
tuations which are the sum of va+iaaa+i, 2wQ 2 a a+1 and 
wAfi(QA,A — Qa+i,a+i)- In this figure and in Fig.|SJ the 
go-representation is used for any value of q in the con- 
vergent region, and Aq (A/i) is numerically determined 
by the configuration dictated method. In the noncon- 
vergent region, the go-representation is fixed at the crit- 
ical point and A/j, is treated as a changeable parameter. 
Effects coming from the higher-order term 0(b) is ne- 
glected. Our numerical results clearly show how nicely 
the analytic condition in Eq. i|38[l holds. Near the criti- 
cal point, one may observe that the difference of s.p. en- 
ergies is almost equal to the quantum fluctuations, i.e., 

£A+x{q)-£A{q) « VA+lAAA+l + 2wQ 2 A A+1 +wA[i(Q A ,A- 

Qa+x,A+x)' Figure 0b) shows the difference of s.p. en- 
ergies €A+i(q) — £a(<z)> the two-body residual interaction 
va+xAAA+x, and the effects from quadrupole deformation 
2wQ 2 AA+1 + wAfi(Q A ,A - Qa+x,a+x), separately. One 
may observe that a competition between one-body po- 
tential ca+i — £a and the two-body residual interaction 
Va+xaaa+x indeed plays a dominant role in determining 
whether or not the concept of CHF mean field is realized. 

Figure |8[a) depicts the ratio h A \ +1 /h A ' > A+1 as a 
function of quadrupole parameter fj, (since the average 
quadrupole moment q has no meaning in the noncon- 
vergent region). One may see that in the convergent 
region the ratio is always positive except for some ex- 
ceptional points, where the deformation fluctuation be- 
tween two-specific orbits Qa,a+i is almost zero as shown 



(l) f c(1) b(1) 1 

h A,A+x = i>A+iM ~ £aM + wA^(Q a +i,a+i - Qa,a)} I \ + > 

With the aid of (|34() . a quantity in 13I|) for the case with n = 1 is expressed as 

h A,A+l _ 1 / 6 A+1 (q ) ~ £a(<7o) - (VA+IAAA+I + 2wQ 2 AA+1 ) 



a (x) b (x) b (x)2 _ a (x)2 I —wA(j,(Qa,a ~ Qa+x,a+x) - 0(b) 
where the higher-order term containing a small parameter &W is written as 

0(6) = -2&« 2 (v A+1AAA+1 + 2wQ 2 AA+1 )+^-(l - 2b {1)2 )Q A:A+1 (Q A+1 . A+1 - Q A , 
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FIG. 7: (a) Energy difference between two specific orbits and 
quantum fluctuations as a function of quadrupole moment; 
(b) s.p. energy difference, two-body residual interaction and 
the sum of deformation fluctuation and quadrupole deforma- 
tion. 



FIG. 8: (a) Ratio of off-diagonal components between the 
first and second iterations for each given /i state. Conv. and 
nonconv. stand for the convergent CHF states and noncon- 
vergent region, (b) Absolute value of off-diagonal quadrupole 
moment as a function of quadrupole parameter fi. 



in Fig. Efb). Near the critical point, the ratio is going 
to reach zero. When one slightly decreases the defor- 
mation A/i from the critical point, the ratio becomes 
negative where the CHF calculation meets a difficulty of 
nonconvergence. After the nonconvergent region, the ra- 
tio becomes positive, and the convergent state appears 
again. From these numerical calculation, it is clear that 
our condition in Eq. I|29l) actually works well in evaluat- 
ing whether the concept of the mean field is realized or 
not in the many-Fermion system. 

We should now comment on the exceptional points dis- 
cussed above. Since the ph component of the quadrupole 
operator IQ^^+il is almost zero, which just corresponds 
to a full alignment state created by the angular momen- 
tum constraint in the cranked HF theory, a new deformed 
state can not be generated by the quadrupole operator 
within the truncated two dimensional subspace. What 
actually happens at these exceptional points attributes 
to the importance of the other part of the s.p. space out- 
side the two dimensional subspace spanned by ipA(qo) 
and ipA+i(qo)- Since our simple analytic understanding 
that uses the 2x2 Hamilton matrix does not work in this 
specific situation, an appearance of exceptional points is 



out of our present discussion. 

V. SUMMARY 

The present work does not include the pairing cor- 
relation. Including the pairing correlation might make 
the situation so far discussed more complicated because 
of many dynamical competitions that exist not only be- 
tween the ph-type two-body residual interaction and the 
HF potential, but also the pp-type two-body residual in- 
teraction and the pairing potential; their cross effects also 
must be considered. Moreover, the two mean fields char- 
acterized by different configurations are mixed up by the 
pairing correlation, and the uv factor introduced by the 
BCS theory obscures the concept of the configuration. 
Although some numerical evidence on an applicability of 
the CHFB theory near level crossing region has been dis- 
cussed in Ref . [T|| , the role of pairing correlation on the 
level-crossing dynamics will be discussed in a separated 
paper. 

In summary, the nonconvergent difficulty in level repul- 
sive region is discussed analytically in the self-consistent 
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mean field theory. It turns out that the CHF mean field 
breaks down when the quantum fluctuations coming from 
two-body residual interaction and quadrupole deforma- 
tion become larger than an energy difference between 
two avoided crossing orbits. Deriving the analytic condi- 
tion, we make it clear that the competition between one- 
body potential and quantum fluctuations mainly coming 
from two-body residual interaction plays an important 
role whether the self-consistent CHF mean field is real- 
ized or not. However, the breakdown of the CHF mean 
field at certain range of the quadrupole deformation q 
does not mean that introducing some quantum mechan- 
ical treatment is necessary, because the constrained op- 
erator Q in the CHF theory is put in by hand. That is, 
the system does not like to be elongated nor contracted 
along a given direction of quadrupole deformation any 
more, but it likes to develop toward a direction chosen 
by itself. Further microscopic investigation is needed to 
answer the very interesting conclusion in Ref. |9j by intro- 



ducing dynamical constrained operators based on the self- 
consistent collective coordinate (SCC) method [32L l33|. 
because a diabolic point related to the level crossing in- 
dicates an existence of missing degree of freedom. 
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